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1 Introduction 

The radiative transport equation (RTE) is a linear Boltzmann equation, which de- 
scribes propagation of energy density in a random medium lIT]|2l. Rigorous anal- 
ysis on the transport limit for the Schrodinger equation with random potential has 
been studied [3 4 5 6 7 8 910111 2|. In particular, the RTE was obtained from 
the Schrodinger equation with time independent [10] and time dependent 111] ran- 
dom potential. In |1T2| . the RTE was derived from the Schrodinger equation with 
random and periodic potential. See also recent review [13] and references therein. 
In the absence of random potential, the semiclassical equations of motion for the 
Schrodinger equation with periodic potential and electromagnetic fields has been 
considered |[Il[T5][l6l[ni[ISllI3- 

In this paper, we consider noninteracting electrons in a periodic potential (Bloch 
electrons) with a weak random potential and apply electromagnetic fields. For this 
system, the RTE is derived in the same way as lfT2l except that we here take elec- 
tromagnetic fields into account. 

This paper is organized as follows. In Sec.|2] a two-scale asymptotic expansion 
is introduced for the Wigner distribution function. In Sec. |3] Bloch functions are 
considered. The Wigner distribution function is decomposed in Sec. |4] and Sec.|5] 
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making use of the Bloch functions. Finally, we derive the RTE and obtain our main 
result Eg. (17911 in Sec.|6l 



2 The Schrodinger Equation 

We consider Bloch electrons in an electric field E and a magnetic field B eR^. 
They are given by vector potential A G and scalar potential (p G /? as 

E = -Va>-^A, B = VxA. (1) 
dt 

The state \lf{t,x) e R (t e R, x e R^) of this system evolves according to the 
Schrodinger equation 

ih^\lf{t,x) = J^Xff{t,x). (2) 
at 

The Hamiltonian is given as 

,^=^ip + eAf + e(p + Utotix), (3) 

where p = —iKVx, and nie are the charge and mass of an electron, and Utot{x) = 
U (x) +V{x). Here U{x) is a c/-dimensional (d < 3) periodic potential and V{x) is 
a random potential. The periodic potential satisfies 

U{z + v) = U{z), (4) 

where v belongs to the lattice L: 

( d 



HjEZ) (5) 



and e 1 , . . . , e^/ form a basis of R"^ with the dual basis e \ e'' defined by e , • e*^ = 
2n5jk- The dual lattice L* is spanned by {e^}. We let C denote the basic period 
cell of L and I^bz denote the Brillouin zone. Note that 

|C|i2Bz = {Inf. (6) 
The random potential has the following properties. 

{V{y)V{y+x))=R[x), {V{q)V{q')) = {2nYR{q')5{q + q'). (7) 
Here (•) denotes ensemble average and the Fourier transform is defined as 

V{q)= f djce-'^-'Y(jc). (8) 
Jr'I 

Because of the {/(I) gauge symmetry, we can set 

9 = 0. (9) 
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We scale the variables as 



t 

t -> -, 
e 



X 

JC — > — , 

e 



where e (> 0) is small. We assume weak electromagnetic fields: 

E-^eE, S eB. 



(10) 



(11) 



Equation ([T]| implies that A is independent of e. Since the cyclotron radius i is 
inversely proportional to we have £ £/e and £ gets much larger than the cell 
size of the lattice L. Furthermore, we assume that the random potential is weak. 
We obtain 



d , , ihe 
1 

+ 



(12) 



where the initial wave function i//£(0,x) G L?{R^) is e-oscillatory. In [12], the RTE 
is derived from this equation in the absence of the vector potential A 
Let us consider the Wigner distribution function associated with i//£. 

W,{t,x,k) = f JL^^'''yrif^(t,x-ey)%{t,x). (13) 

jRd [27tf 

Whenever necessary, we use the fact that G is uniquely decomposed as 

k = q + ii, qeQBZ, pi el*. (14) 
Note that Ws{t,x,k) and its symmetric version 

[Inf 

have the same weak limit as £ — !■ [20J. We obtain the time evolution of as 



—We{t,x,k) + 



hk „ me e / , s „ , „ A 



ihe 



l£ 



— (Ait,x) + -A{t,V,))-A{t,V,) 



Wsit,x,k) 



J_ £ e'^''-/-u{p')[Ws{t,x,k-p')-Wsit,x,k)] 



dk' 



\W^{t^X,k-k!)-We{t,X,k)\ , (16) 



where we used 



\lfe{t,x) — \if^{t,x) 0, {i=\,...,d), 



(17) 
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and introduced 



U{y) = I e'>->'C7(Ai), {/(M) = ^ / dye-'>->'[/(.y) 
,i7l* \C\ Jc 



Note that we have 

Let us define 
Noting that 

we obtain 
d 



1 



h 



H{x,lJ.,q) = 



2mp 



n 



dA{x) d 



= AC(0,. .. ,0,^,0,. ..,0)), 



+ [Vi,H{x,ll,q) ■ V, - V,H{x,ll,q) ■ V^] = 



We note that the term —VxH is responsible for the Lorentz force: 



-V,(^-A) + 0(|A|2), 



-— V,(fc-A) = -^p.V,)A + fcx(V,xA)] 
me nie 

n 

where hk = m^v + eA and v — dx /df . The left-hand side of Eq. 
pressed as 



(18) 
(19) 
(20) 
(21) 



(22) 



d_ 
dt' 



1, 



1, 



1. 



1, 



This expression was first obtained by Kubo II21L 
We introduce a two-scale expansion for W^: 

W^it ,x,k) = WQ{t ,x,z,k) + ^/eWy {t,x,z,k) + eW2{t ,x,z,k) 



(23) 
can be ex- 

Wi;{t,X,k). 

(24) 
•••• (25) 



We assume that Wq is deterministic and periodic with respect to z = x/e. We re- 
place 

V,^V, + -^V,. (26) 
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By collecting tenns of 0(e ^ ), we have 

where the skew symmetric operator ££ is given by 



^f{z,k) = V^H{x,^,q)-VJ{z,k) + —AJ{z,k) 

-4 I e^'-'u{^')[f{z,k-^')-f{z,k)] 



hk e , / X 

\ A{t,x) 

me me 



ih 

■^zf{z,k) + —AJ{z,k) 
2me 



1 ^ e''''-i/(M')[/(z,^-M')-/(^,^)]- 



By collecting terms of 0{e ^/^), we have 



By collecting terms of 0(1), we have 

d_ 
dt 



^Wo + (V^// • V;, - V^H ■ V^) Wo = -^W2 - ^V, ■ V,Wo{t,x,k) 



+ ^ 7^e*'-^y(/:') [Wi{t,x,k-k')-Wi{t,x,k)] . (30) 



(27) 



(28) 



(29) 



{2n) 



3 The Bloch Functions 

To obtain the eigenvalues and eigenfunctions of the operator _Sf , we consider the 
following eigenproblem. 



^0V, + M(.,x)) +Uiz) 



'i>ma{z,q)^En{q)<Pna{z,q)- (31) 



Here, a labels degenerate energy levels of E„ with multiplicity r^: a = 1, . . . , r^. 
We assume that there is no level crossing and 

Ei{q)<E2{q)<---. (32) 

The parameter q&R'^ labels eigenvalues of the translational operator T {— e^'^s 
veL): 

Tcl^(z.q) = «>(z + V.ty) = ^{z,q). (33) 

Note that q moves inside the Brillouin zone (I^bz)- It is, however, convenient to 
extend ^{z,q) to in q with L* -periodic. The eigenfunctions ^ma{z,q) form a 
complete orthonormal basis in L^{C): 

f dz 

(*ma, = 1^ *ma(z, (z, = S,njSap ■ (34) 
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We rewrite the eigenproblem using the periodic function 

^iz,q)=C~"^-'<P{z,q). (35) 

We obtain 

2me 



^ (^V-^+hq + eA{t,x) ] +U{z) 



<l>ma{z,q)=Eni{q)^ma{z,q). (36) 



By differentiating the equation with respect to qj, we obtain 



Thus we have 

dE,n „ „ .h^ ( d'Pina ^ \ , 

0„„,0„R = -I ,<P„R H J 



dqj ^ \ dzj J ' 



Similarly, we also have 



' J " ~iq^ ' ^'li^^^^' {^l^ma , *„/3 ) , (39) 

where 

/" dz 

(A, , *„/j ) = 1^ A, (f , z) <P,„a (z, ^) *„j3 (z, ^) . (40) 
The Bloch functions satisfy the following orthogonality relations II12I . 

77— £ / dq(p„,aix,q) 4>ma iy,q) = S (v - x) , (41) 
J^^dx(pjaix,q)i>,„p{x,q') = 5j„,5ap5{q - q). (42) 

4 Decomposition of Wq 

Let us define the z-periodic functions Q,„„{z,lJ.,q), jJ. eL*,q e Qbz by 

Q"l (z, At , ^) = ^ ^ e'(^+'') >*„„ (z - y, q) <P„p (z, ^) . (43) 
These functions satisfy the following orthogonality relation. 

I /^ef (z,M,^)e^''''(z,M,^) = 5y,„5„„,5p/5'- (44) 
They are eigenfunctions of 

J^Q^P {z,^l,q) = ^ [E,„ {q) - E„ {q)] Q^i (z, M , ^) , (45) 



(38) 
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with k = q + ix. 

Let us write Qmm as Q,"'^ . Equation ( 145 b implies that ker^ is spanned by Q'^ . 
By Eq. ( 127b . Wo(f,Ji:,z,^) may be decomposed as 

W'o(f,Jf,z,^) = W'o(f,x,z,^ + ju)= £ {M,„(f,x,^)}„j3G,^'^(z,jU,^), (46) 

m.a.fi 

where ^) is a r,„ X r,„ matrix. 



5 Decomposition of Wi 



Let us look at Eq. ( [29] l. In general, Wi is not periodic in the fast variable z- Hence, 

is fi 

dv 



instead of Qmn , we use the basis functions 



P^i {z, Ai , ^, ^o) = ^ 1^ e'(^+^') >'*,„„ (z - y, q) <S>„p (z, ^ + ^o) , (47) 

where z E R"^ and q,qo G i^BZ- The functions P"/^ are quasi-periodic in z: 

C^z + v,M,^,^o) =Cf (z,M,^,?o)e-'"-^°, (48) 

where V € L. Using Eq. (fT9] l, Eq. (|34] |. and Eq. (|42] |. we obtain the following 
orthogonality relation. 

I ■^P„lPiz,H,q,q')Pfi^'{z,H,q,q") = 5„,j5Aa'5pp,5{q' - q"). (49) 

The functions P"n are also eigenfunctions of Jff. 

^P^i{z,^,q,qo) ='-[E,n{q)~En{q + qo)]P^!{z,lX,q,qo). (50) 



We write W\ in this basis as 

Aq' 

mnap 



Wi{t,x,z,q + ^)= £ / ^ri^P{t,x,q,q')P,f,P{z,^,q,q'), (51) 



where z e 7?'', q e X2bz, and G L* . We plug this into Eq. ( |29] |, multiply P"'^ (z, pL,q,qo), 
sum over jx £ L*, and integrate over z&R'' ■ The Ihs becomes 

£ I AzPffiz.liA^qo)^^ £ / ^77«/(f,X,^,^')C/(2,M,^,^') 

=^ [^^./l^) -■E/(^ + ^o)]'7,7'^(f,-x,^,?o)- (52) 

We define ifTIl 
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By taking the sum over jj., the rhs becomes 



(27r)(^-l)/2r««'(^,^ - q' - n'){u„,{t,x,q- q')}a'P'i>n.p'{z,q- q') 



SjmSaa'{Um{t,X,q)}a'l}'e''^'l''^^'^''-^„,B'(z,q) 



*Zj3(z,^ + ^o) 



1 V 



y(-^o-M)L7;""'(^,^+^o + M){"Kf>-^,^+^o)} 



a'/3 



(54) 



Therefore we obtain 

aPf ^ ^ ^BZ V- V{-qo-ll) 

r,., [t,x,q,qo) ^2n)i''+m L.^^EM + qo)-EM) + i^ 

X '^{ui{t ,x,q + qQ)}a'pT^,"' {q,q + qQ + ll) 
a' 

r dzdk' i,,'. yik')Y.p'{uj{t,x,q)}ap,(Pip {z,q + qo)^jpi{z,q) 
jR2^{2nY^ Ei{q + qQ)-Ej{q)+i^ 

where t, (> 0) is infinitesimally small. 



(55) 



6 Time Evolution of Wq 

Next let us look at Eq. OOl) . We multiply (z, /i , integrate both sides over z, 
and sum over /J . 



|„-„+A(*+l.(,,,)).v,»'„-^£(*, + l.,(.. 



— ^ A, f,v^ pyo 

me c'X; We ' 



me 
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dk' 



(56) 



The first integral on the rhs vanishes is skew symmetric and Qj ^ G ker^). 
6.1 LHS 

The Ihs of Eq. ( l56l) is calculated as follows. The first term is easy. 



=-^^{Uj{t,X,q)}al}. 



(57) 



The second term is calculated as follows. 

I./,^erfc.,,)(*, + |AK,))(^-iA,(v, 



X 



-4 



= — ;^ Oj„,5^p, -I I — ^ — ,^ja 1 +TAz(ji:)5a„/ 



5x/ 



Therefore, the Ihs of Eq. (|56] | is written as 

5 1 e 

— {Uj{t,x,q)}ap + jVqEj{q) ■ V ^{uj{t ,x,q)} ap + - 

^ L (^V,^-^y^(z)*ya'(z),'^;a(z)) {uj{t ,x,q)}a'p 



(58) 



(59) 
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6.2 RHS 

Let us calculate the rhs of Eq. ( 156b . After taking ensemble average, the rhs is given 
by the sum of two matrices /i and h where 

X {V{l<')Wi{t,x,z,k-k')), (60) 
X {V{k!)Wx{t,x,z,k)). (61) 



Using Eq. (ISTT i and Eq. (|55] |. /i is written as 

/i=/ii-/i2, (62) 

where 

(2;r){rf+l)/2 ^,1-^^ En iq-q' + q") - E„, {q-q') + 

X £{M„(f ^ - ^' + q")}a"P'T„i""iq -q',q-q'+q" + m") 
a" 

xP^'/{z,^^^',q-q',q"), (63) 

and 

>< l,,^^''''"iV(k')V{kO)P,i^'iz,^-^\q-q\q'^) 

Ll}"{umit,x,q- q')}a'p"'Pnp'{z',q - q' + q")^,„p" iz',q - q') 

^ En{q-q' + q")-Em{q-q')+it, 

(64) 



We obtain 



X + + (65) 
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Similarly we also obtain lu- 



R{-k') 



ih jRd 2n ^ Ej {q) - £„, {q-q') + it, 

X Y,Mt,x,q-q')}a'p''Tp'iq,q-k')TP'jl^{q-k',q), (66) 
a'P" 



where we used 



(67) 



We have 



dq' R{-q'-ll') 



(q) - Emiq ~ q') + it, 



£ Tjf,f{q,q- 4 - \i!)T<"-"{q - ^' - [i! ,q){uj{t ,x,q)\a"p 



L '^fm 4 - l^'){umit,x,q- q)}a'fi"T!;,j{q -q'-ll',q) 

a'P" 



(68) 



In the same way, we have 



{h}ap = l^L L 



dq' R{-q'-lx') 



'^^ n pj^L* -'^^ 2^ ^« i^ + l')- E}{q) + it, 

£ T^f {q,q + q' + ^l'){u„{t,x,q + }a"/J' T^j^ {q + q' + ^l\q) 
a"P' 



£ {uj {t,x, q)}ap"Tf^^' {q, q + 4 + ju') T,'^-'^ {q + q' + ^l\q) 
P'P" 



(69) 
(70) 



Thus we see the relation 

h=n, 

where * denotes the Hermitian conjugate. By summing /i and I2, and using the 
relation 

(71) 



T!;'{q,q') = T;Piq\q), 



the rhs of Eq. ( 156b becomes 



Tjm {q,q - 4 - ll')um{t ,x,q - q')T,„j{q - q - ll\q) 
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+ 



Ej{q) -Emiq-q')- Ej{q) -Em{q-q')+ 
Tjm{q,q- 4 - IJ-')T*„{q, q-q' - li')uj{t,x,q) 

Ej{q)-E,„{q-q') + i^ 
uj{t,x,q)Tj^{q,q- q' - n')T*^{q,q- q' - jx') 



Ej{q)-En[q-q')-i^ 
We change the variable asq' —>q — q' and obtain 

2mTj,n {q,q' - ^l')um{t ,x,q')T„j{q' - ^l' ,q)S {E j{q) - E,„{4)) 

Tj^jq.q' - lj')TjJq.c/ - Ij')llj(t.x.q) 
^ Ej(q)-E„{q') + i^ 

_ Uj{t,x,q)Tjm{q,q' - ^l')T*Jq,q' - 
Ejiq)-E^{q')-i^ 



(72) 



(73) 



6.3 Radiative Transport; Equation 

We define a vector Vj{q), a matrix Mj{q), and a superoperator /il as 



{Mi{q)}ap = y-Vj{q)-A4>jp,4>ja j , 
Al : Uj{q) ^ [Mj{q),uj{q)] ^ Mj{q)uj{q) - uj{q)Mj{q). 

Furthermore let us define the following superoperators. 



A{q,q'):Uj{qo)^ ^ R{-q' - ^')5 {Ej{q) - E^)) 

Tjm{q,Ct - ^i'')Um{qo)Tmj{ct - ^l'' a), 
-1^ ^ /• d^' 



in ii'sL- 

Tj,„{q,q' - ^')T*Jq,q' - ^l')uj{q) 
Ejiq)~E,„{q')+i^ 

Uj{qo)Tj,n{q,q' ~ ^l')T*Jq,q' - mQ 
Ej{q)-E^{q')-i^ 



(74) 

(75) 
(76) 



(77) 



(78) 
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Finally,by Eqs. OOl, (ED, and (|73]l, we obtain the RTE as 



—Uj{t,x,q)+Vj(q)-VxUjit,x,q) + + fls)uj{t,x,q) 

= [ Aq'A{q,q')uj{t,x,q'). (79) 
JBZ 

Note that the term for /1l stems from the Lorentz force. The terms for /is and A 
stem from the random potential. 
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